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ABSTRACT 

Both analytical and numerical works show that magnetic reconnection must occur in hot accretion 
flows. This process will effectively heat and accelerate electrons. In this paper we use the numerical 
hybrid simulation of magnetic reconnection plus test-electron method to investigate the electron accel- 
eration and heating due to magnetic reconnection in hot accretion flows. We consider fiducial values 
of density, temperature, and magnetic parameter (3 e (defined as the ratio of the electron pressure to 
the magnetic pressure) of the accretion flow as no ~ 1 6 cm~ 3 , Tj? ~ 2 x I0 9 K, and (3 e = 1. We find 
that electrons are heated to a higher temperature T e — 5 x I0 9 K, and a fraction r\ ~ 8% of electrons 
are accelerated into a broken power-law distribution, dN^) oc 7~ p , with p w 1.5 and 4 below and 
above ~ I MeV, respectively. We also investigate the effect of varying (3 and no. We find that when 
P e is smaller or no is larger, i.e, the magnetic field is stronger, T e , n, and p all become larger. 
Subject headings: acceleration of particles - accretion, accretion disks - black hole physics - magnetic 
fields 



1. INTRODUCTION 

Magnetic reconnection serves as a highly efficient en- 
gine which converts magnetic energy into plasma thermal 
and kinetic energy. The first theoretical model for mag- 
netic reconne ction was proposed by Sweet and Parker 
(lParkerlll958fl . Ever since then, theoretical studies along 
this line have been reported in various papers includ- 
ing analytical theories and numerical simu lations (jParkerl 
[1961 iDrake et al.ll2005l : lLitvinenkd[T99l . 

Magnetic reconnection has been widely used to explain 
explosive phenomen a in space and laborato ry plasmas 
such as solar flares (Somo v fe KosugH I1997T), the heat- 
ing of solar corona (jCargill fc KlimchukllI997D. and sub- 
storms in th e Eart h's magnetosphere (jBieber et al.llI982 : 
I Baker et al.l 1200 25. Magnetic reconnection also plays 
an important role in high-energy astrophysical environ- 
men ts such as the magnet ized loop of the Galactic cen- 
ter dHevvaerts et al.l If 9881). iets in active gal a ctic nuclei 
(Scho pper et alJ 119981 : lLarrabee et alJ 120031 : iLvutikovl 
120031 ). and t rav bursts ([Michel 1119941 ). For example, by 
assuming that particle acceleration is due to direct cur- 
rent (DC) electric field and is balanced by synchrotron 
radiative losses, Lyutikov (2003) found that the maxi- 
mum energy of electrons is ~ 100 TeV. Larrabee cal- 
culated the self-consistent current density distribution 
in the regime of collisionless reconnection at an X-type 
magnetic neutral point in relativistic electron-positron 
plasma and studied the acceleration of non-relativistic 
and relativistic electron populations in three-dimensional 
tearing configurations. 

But so far there is very few work on the magnetic 
reconnection in hot accretion flows. In accretion flows, 
the gravitational energy is converted into turbulent 
and magnetic energy due to the magnetorotational 
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120081 ) and three-dimensional magneto hydrodynamic sim- 
ulations of black hole accretion flow faawlev fc Balbusl 
12001 iMachida fc Matsumotol 12001 iHirose et all 120041) 
have shown that magnetic reconnection is unavoidable 
and plays an important role in converting the magnetic 
energy into the thermal energy of particles. In fact, 
turbulent dissipation and magnetic reconnection are 
believed to be the two main mechanisms of heating and 
accelerating particles in hot accretion flows. In the case 
of large (3 ((3 is defined as the ratio of the gas pressure 
to the magnetic pressure), the magnetic reconnection 
is even likely to the domi nant mechanism of electrons 
heating and acceleration (IQuataert fc Gruzinovl 119991 : 
iBisnovatyi-Kogan fc Lovelace! Il998f ). On the obser- 
vational side, flares are often observed in black hole 
systems. One example is the supermassive black hole 
in our Galactic center, Sgr A*, wher e infrared and 
X-ray flares occur almost ever y day (e.g., Baganof f et alJ 
I200lt iHornstein etaD 120071 : IDodds-Eden et all 12009). 
To explain the origin of these flares, we often require 
that some electrons are transiently accelerated into a 
relativistic power-law distribution by magnetic recon- 
necti o n within the accretion flow (e.g., B aganof fet al.l 
pOOll lYuan et all l200l l200l IDodds-Eden et alJ 120091 : 
lEckart et al.l 12003 ) or in the corona above the accretion 
flow (e.g., Yuan et al. 2009), and their synchrotron or 
synchrotron-inverse-Compton emission will produce the 
flares. 

Regarding the study of the electron acceleration by 
reconnection, recently many works have been done on 
the relativistic magnetic reconnectio n with the approach 
of Par ticle- in-Cell (PIC) simulation. IZenitani fc Hoshind 
(2007) presented the development of relativistic magnetic 
reconnection, whose outflow speed is on the order of the 
light speed. It was demonstrated that particles were 
strongly accelerated in and around the reconnection re- 
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gion and that most of the magnetic energy was converted 
into a "non-thermal" part of plasma kinetic energy. Mag- 
netic reconnection during collisionless, stressed, X-point 
collapse was studied usi ng kinetic, 2.5D, fully rela tivis- 
tic PIC numerical code in lTsiklauri fc Haruki (200?j) and 
they also found high energy electron s. PIC electrom ag- 
netic relativistic code was also used in lKarlickyj (|2008t ) to 
study the acceleration of electrons and positrons. They 
considered a model with two current sheets and periodic 
boundary conditions. The electrons and positrons are 
very effectively accelerated during the tearing and coa- 
lescence processes of the reconnection. All these work 
investigate the circumstance of the Sun. 

In this paper we use the numerical hybrid simulation 
of magnetic reconnection plus test-electron method to in- 
vestigate the electron heating and acceleration by mag- 
netic reconnection in hot accretion flows. We would like 
to note that the results should also be applied to corona 
of accretion flow and jet, if the properties of the plasma in 
those cases are similar to what we will investigate in the 
present paper. Self-consistent electromagnetic fields are 
obtained from the hybrid code, in which ions are treated 
as discrete particles and electrons are treated as mass- 
less fluid. Thereafter, test electrons are placed into the 
fields to study their acceleration. The input parameters 
include density, temperature and magnetic field of the 
accretion flow. The values of density and temperature 
of astrophysical accretion flows can differ by several or- 
ders of magnitude among various objects. H ere we take 
the a ccretion flow in Sgr A* as reference (see I Yuan et al.l 
120031 for details), but we also investigate the effect of 
varying parameters. For the strength of the magnetic 
field, we consider mainly (3 e = 1. Here [3 e is defined as 
the ratio of the electron pressure to magnetic pressure, 
(3 e = Pe/Pmag- Note that what we usually use in hot ac- 
cretion flows is the ratio of the gas pressure and the mag- 
netic pressure, /? = -P ga s/-fmag- The hot accretion flow 
is believed to be two-temperature , with Tj ;g> T e . F or a 
reasonable value of T t = 10T e (ref. lYuan et al.ll2003D . the 
above /3 e corresponds to the usual j3 = 10. This is the 
typical value in the hot accre tion flow as shown by MHD 
numerical simulations (e.g., iHirose et al.l 12004 ). Given 
that the value of (3 is ver y inhomogenous in the accretion 
flow (jHirose et al.l 120041) . we also consider [3 e — 10,0.1, 
and 0.01. 

In Section 2, we use the numerical hybrid simulation 
of magnetic reconnection to get the structure of self- 
consistent electric and magnetic fields. In Section 3 we 
use the obtained electromagnetic fields as the background 
of test electrons to investigate the electron acceleration. 
Finally we summarize our results and discuss the appli- 
cation in interpreting the flares in Sgr A* in Section 4. 



2. THE NUMERICAL HYBRID SIMULATION OF 
MAGNETIC RECONNECTION 

We use the 2.5-D (two-dimensional and three compo- 
nents) hybrid simulation code of Swift (jSwift Ill995lll996f ) 
to model the driven steady Petschek mode magnetic re- 
connection, except that we use a rectangle coordinate 
system. Quasi charge neutrality is assumed, which can 
be written as rii « n e — n. 

The momentum equation for each ion as a discrete par- 
ticle is given by 



^ = (E + vxB)-!/(ui-u e ), (1) 

where v is ion velocity, E is the electric field in units of 
ion acceleration, B is the magnetic field in units of the ion 
gyrofrequency, and v is the collision frequency. A finite 
numerical resistivity corresponding to the collision fre- 
quency v is imposed in the simulation domain (X point) 
to trigger the reconnection. u e and are the bulk flow 
velocities of electrons and ions. In particular, is cal- 
culated by CIC (coordinate in cell) method from v and 
the electron flow speed u e is calculated from Ampere's 
law: 

u e = Ui - (y x B/nu e). 

The electric field is derived from electron momentum 
equation 

„ , 1 m e du e 

E = — U e X B — v(\L e — Uj) V-Pe TT- (2) 

ne mi at 

where P e is obtained from the adiabatic equation 
P e (m e n) 5 / 3 = constant, n is the density, and e is the 
charge of electron. In equation (2), m e <C rm, so the last 
term can be omitted, and thus we can define electric field 
E with equation (2). 
We use Faraday's law to update the magnetic field 



dB 

~dt 



V xE. 



(3) 



The velocity of ion is updated at half time step using 
a leapfrog scheme with a second order accuracy. The 
magnetic field and the particle positions are advanced 
at full time step using an explicit leapfrog trapezoidal 
scheme. 

Initially a current sheet separates two lobes with anti- 
parallel magnetic field in the x-direction. The current 
sheet is located along z=0.5L z . The sum of thermal and 
magnetic pressure is uniform across the initial current 
sheet. Initial conditions are: 



Bo = Bq x t&nh.(z/5)e x 



no = «oo + 
T e = n T . 



re cosh 2 (z/ (5) ' 



(4) 



Here 8 is the half width of the initial current sheet. Bq x , 
rioo is the initial physical variables near the boundaries 
z—0, L z . The initial electron temperature To is constant, 
re is the ratio of the boundary density to current sheet 
density and we set re = 4. 

At the boundaries z~0 and L z , B z is set to be zero, 
and there is a small inflow velocity u x on the two bound- 
aries. The inflow is continuously imposed on the bound- 
aries throughout our simulation. Spatially, the inflow is 
uniform along the boundary. 

For the parameters of the black hole accretion flow, we 
take the accretion flow in Sgr A* , the superm assive black 
hole in our Galactic center, as an example (|Yuan et al.l 
2003). The characteristic density, temperature, and f3 e (= 
Pe/P mag ) are n e ~ 10 6 cm- 3 ,T e - 10 9 3 K, and (3 e = 1.0, 
respectively. Given the inhomogeneity of the accretion 
flow we also consider (3 e = 10, 0.1, and 0.01. 
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In the simulation the cell size in ^-direction is chosen to 
be Az = Ao (Ao = c/uj p io is the ion inertial length), and 
the cell size in the ^-direction is Aa; = 2Az. The length 
of the simulation region is L x = lOOAcc and L z = lOOAz. 
The magnetic field B is in units of ion gyrofrequency 
ilo, time is in units of ^q 1 , and v and u are both in 
units of Va (Va is the characteristic Alfvenic velocity on 
boundaries z=0 or L z ). 

With the above units, we have chosen 5=0.45 and 
n. ce «=150 particles per cell. The spatial profile of the 
resistivity imposed in the simulation corresponds to a col- 
lision frequency v = ^ocxpj— [(x — xq) 2 + (z — zo) 2 ]/i5 2 }. 
Here Xq and zq are the middle point of simulation do- 
main, with xq—0.5 L Xl zo=0.5 L z , and vq =0.02. The 
time step to advance the ion velocity is chosen to be 
At = O.IQq 1 . The evolution of the configuration of the 
magnetic field is shown in Figure 1. Induced electric field 
E y (which is vertical to x-z plane) plays an important role 
in electron acceleration. Spatially, the inductive electric 
field appears mainly in the two pile-up regions (on the left 
and right hand sides of the X-point) on the current sheet, 
having the shape of two circular spots of high value. As 
the reconnection goes on, the outflow together with the 
magnetic field moves out from the center of the simu- 
lation box because of the magnetic frozen-in condition. 
Along with it, two spots of inductive electric field also 
move out slowly along the current sheet in opposite di- 
rections. Temporally, it varies little during the electron 
acceleration and can be regarded as steady. Figure 2 
shows the evolution of the induced electric field. 



*/L AT TIME 1 CO 



Fig. 1. — The evolution of configurations of the magnetic field 
during the reconnection process. 



3. ELECTRONS HEATING AND ACCELERATION IN HOT 
ACCRETION FLOWS 

Magnetic reconnection is of course intrinsically a tran- 
sient process. But in our paper we treat the electron 
heating/acceleration in reconnection as approximately 
steady. Specifically, from Figs. 1 & 2 we see that 
the electromagnetic field almost remains unchanged from 
t=75 to 125 and reaches a steady state. We investigate 
the heating/acceleration of electrons during this interval. 
This approximation is reasonable because the timescale 
of electron heating/acceleration is very short compared 
to the timescale of the reconnection process. 



x/L AT TIME 75 



AT TIME 100 




Fig. 2. — The evolution of configuration of the inductive electric 
field during the reconnection process. 

Electron dynamics in magnetic reconnection has been 
previously studied using analytical theories and test 
particle calculations in different magnetic- and electric 
field configurations. (|Litvinenkolll996t iDeeg et al1 | 1991 
Browning k Veksteinl 120011 : iWood fc NeukirchT 12005 : 
Birn &: Hessd 19941 ) . In such studies, the electron orbits 
are calculated in given electromagnetic fields. Electrons 
are m ainly accelerated by the v x B reconnection electric 
field. iLiu et all (|2009i ) performed test-particle simula- 
tions of electron acceleration in a reconnecting magnetic 
field in the context of solar flares, and found that the 
accelerated electrons present a truncated power-law dis- 
tribution with an exponential tail at high energies, which 
is analogous to the case of diffusive shock acceleration. 

We use the test-electron method to investigate the 
acceleration. The electromagnetic fields obtained from 
2.5D hybrid simulation of magnetic reconnection is used 
as the background. One advantage of this method is 
that it saves a lot of computation time and enables us to 
perform more calculations. 

It is hard to study the electrons movement in the time 
scale fly 1 and space scale Ao in hybrid simulation. So 
the cubic spline interpolation is used to get electric and 
magnetic fields with higher time resolution ( 0.001 FIq 1 , 
one step is divided into 100 steps, the step in hybrid 
simulation is 0.1 ^q" 1 ), suitable for the calculation for 
electron dynamics. And we also use the values of electric 
and magnetic fields on four grids around an electron to 
calculate the local electric and magnetic fields where the 
electron locates. The difference method has the following 
form. 



F(x ) z) = fu0--dx)(l-dz) + M(l 
+/i2<ix(l — dz) + j 2 idxdz. 



dx)dz 



Here the function F can be any component of the electric 
and magnetic fields; /n 21,12,22 are their values on four 
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grids around the electron; dx and dz are the distance 
from the electron to the nearest grid in x and z direc- 
tion, respectively. Such a treatment will rule out arti- 
ficial high frequency disturbance. The accuracy of such 
difference method is of second order and also satisfies the 
requirement of our simulation. The electron movement 
should satisfy the Newton equation (the relativistic effect 
is considered): 



^=e[E+(p/ 7 )xB]. 



(5) 



Where p is the electron momentum and 7 is the Lorentz 
factor. 
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Fig. 3. — The electrons energy distribution before and after ac- 
celeration in the case of f3 e = 1.0. The dotted line denotes the 
initial (Maxwellian) distribution. The solid line shows the distri- 
bution after acceleration. It can be fitted by a new Maxwellian 
distribution (shown by the dashed line) plus a hard tail. 
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Fig. 4. — The distribution of accelerated electrons ("hard tail" in 
Fig. 3) when f3 = 1.0. It is fitted by a broken power-law. 
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The size of the simulation box is [— 50Ao, 50A ] x 
[— 25Ao, 25Ao]. Because of symmetry and in order to 
reduce the calculation, test electrons are added is in 
[0, 50Ao] x [0, 25Ao] the region only. The number of the 
test electrons is 65536. Their initial velocity distribution 
is Maxwellian and the spacial distribution is homoge- 
neous. The motion of test electrons is controlled by the 
electric and the magnetic fields obtained in §2. 

The resulted new distribution after the acceleration in 
the case of (3 e — 1.0 is shown in Figures 3 & 4. We 
see that the distribution can be described by a new 
Maxwellian one with higher temperature plus a hard 
tail which can be fitted by a broken power-law, i.e., the 
electrons are heated and accelerated after the magnetic 
reconnection. The temperature of heated electrons in- 
creases from 10 9 ' 30 -ftT to 10 9 65 K. The broken power-law 
is described by dN(E) oc E~ p dE, with p = 1.47 between 
1 MeV and 10 MeV and p = 4 between 10 and 40 MeV. 
Most electrons are heated rather than accelerated, and 
about 7.94% of the total test electrons are accelerated 
above lMeV and 0.93% are accelerated above lOMeV. 

We have also calculated the cases of other values of 
(3 e , density n e , and initial electron temperature T®. The 
results are shown in Figures 5 & 6 and Table 1. We can 
see that the results are qualitatively similar. Electrons 
are heated to a higher temperature, and some of them 
are accelerated. The accelerated electrons can be well 
described by a broken power-law distribution. When 



Fig. 5. — The electrons energy distribution before (dot-dot-dot- 
dashed line) and after acceleration in the cases of fi = 10.0 (dashed 
line), /3 = 0.1 (dotted line), and /3 = 0.01 (solid line). 

the magnetic field becomes stronger, i.e., (3 e becomes 
smaller, electrons are more effectively heated and accel- 
erated, namely the "new" temperature becomes higher, 
and more electrons are accelerated. This implies that 
more magnetic energy is converted into the kinetic en- 
ergy of electrons with the increase of the magnetic field. 
When /3=0.1, some electrons can even be accelerated to 
more than 5000 MeV. When /3=0.01, however, the mag- 
netic field B is too strong. In this case, the accelerated 
electrons reaching a certain energy will be deflected by 
the strong magnetic field and run out of the acceleration 
region, thus further acceleration is restrained. Although 
more electrons can be heated to higher temperature and 
accelerated to higher energy when (3 decreases from 0.1 
to 0.01, the highest energy that the electrons can reach 
when (3 = 0.01 is smaller than that when (3=0.1. Another 
interesting result is that the "low-energy" part of the bro- 
ken power-law, where most accelerated electrons locate, 
becomes softer with decreasing (3 e , with the spectral in- 
dex p ~ 0.63,1.47,1.46 and 1.91 for e = 10.0,1.0,0.1, 
and 0.01, respectively. 

When the density n e increases (but (3 e and T° remain 
the same), the effect is similar to the case of decreas- 
ing (3 e , namely the electrons are more effectively heated 
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TABLE 1 

Results of electrons acceleration for different initial temperature, density and fie 



r»(K) 


n e (cm 3 ) 




Te(K) 


power-law form 


acceleration ratio 


10 9.3 


10 6 


10 


10 9.55 


dN(E) oc E-° G3 (1 -14 MeV) 
rIN(F\ W.~ 7.59 /ivi oo y„vi 

til V I _C/ I (_X J_y I l^i-ZO ivlc V ) 


2.24%(E>1 MeV) 
u. ov /o iu ivie v i 


10 9.3 


10 6 


0.1 


10 9.85 


dJV(B) oc B- 1 - 46 (l-17 MeV) 
dN(E) oc £ _5 - 05 (17-46 MeV) 


44%(E>1 MeV) 
3.2%(E>10 MeV) 
0.5%(E>1000 MeV) 


10 9.3 


10 6 


0.01 


10 10.5 


£ZiV(£) ex E- 1 - 91 (5-20 MeV) 
dN{E) oc £- 3 - 80 (20-100 MeV) 


63.1%(E>5 MeV) 
28.8%(E>10 MeV) 


10 9.3 


10 6 


1 


10 9.65 


dN(E) oc S- 1 - 47 (l -10 MeV) 
dN(E) oc £" 4 - 03 (10-40 MeV) 


7.94%(E>1 MeV) 
0.93%(E>10 MeV) 


10 9.3 


10? 


1 


10 9.75 


dN(E) oc E- 1 - 83 (l-10 MeV) 
dN(E) oc £- 3 - 25 (10-35 MeV) 


17.3%(E>1 MeV) 
1.1%(E>10 MeV) 


10 B.8 


10 b 


1 


10 9.2b 


dAf(E) oc S- i Ba (0.5-2.1 MeV) 
dN(E) oc £T 4 - 84 (2. 1-4.5 MeV) 


5.0%(E>0.5 MeV) 
1.6%(E>1 MeV) 
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Fig. 6. — The distribution of accelerated electrons in the cases 
of /3 = 10.0 (solid line), /3 = 0.1 (dotted line), and (3 = 0.01 (dash- 
dotted line). All of them can be fitted by a broken power-law 
(dashed line). 

and accelerated. This is because larger density implies 
stronger magnetic field which is the dominant factor de- 
termining the heating and acceleration. When the initial 
temperature Tj? decreases (from 10 9 3 K to 10 8 5 K) but 
density and f3 e remain unchanged, we find that the heat- 
ing and acceleration become less efficient. The "heated" 
temperature is lower, and less electrons are accelerated. 
Only a fraction ~ 5% of electrons can be accelerated to 
more than 0.5 MeV. This is because a lower T° corre- 
sponds to a weaker magnetic field. But surprisingly, the 
low-energy part of the power-law distribution becomes 
softer which seems to indicate that magnetic field is not 
the only factor determining the acceleration. 

4. SUMMARY AND DISCUSSION 

In this paper we use hybrid simulation code of mag- 
netic reconnection plus test-electron method to study the 
electrons heating and acceleration by magnetic reconnec- 
tion process in hot accretion flows. The self-consistent 
electromagnetic fields are obtained from the hybrid sim- 
ulation and they are then used in the test-electron cal- 
culation. The fiducial parameters of the background ac- 
cretion flows we adopt are density n e = 10 6 cm -3 and 



electron temperature T e ° — 2 x 10 9 3 K (ions temper- 
ature Xi = 10T"). These values are taken from the 
accretion flow in Sgr A*, t he supermassive black hole 
in our Galactic center (ref. lYuan et al.l [20031 . For the 
strength of the magnetic field in the accretion flow, we 
set P e (= Pe/Pmag) = 1-0. This is the typical value ac- 
cording to the MHD numerical simulation of hot accre- 
tion flow. We find that the electrons are heated and 
accelerated by the reconnection. The new distribution 
can be fitted by a Maxwellian distribution plus a broken 
power-law. The temperature is increased to 10 9 65 K and 
the power-law is described by dN(E) oc E~ 1A7 between 
1 and 10 MeV and dN(E) oc E~ 4 between 10 MeV and 
40 MeV. The fraction of electrons with energy above 1 
MeV is ~ 8%. Given that the magnetic field and density 
in accretion flow are inhomogeneous, we also consider 
e = 10,0.1, and 0.01 and n e = 10 7 c m - 3 . We find that 
the results are qualitatively similar, namely the distribu- 
tion of electrons can be fitted by a Maxwellian one with 
a higher "heated" temperature plus a broken power-law. 
When the magnetic field is stronger (/3 e is smaller or n e 
is larger), the heating and acceleration become more ef- 
ficient. The "heated" temperature is higher and more 
electrons are accelerated. In addition, the "low-energy" 
part of the broken power-law becomes softer (refer to 
Table 1). 

In our work we assume ideal MHD for the hybrid sim- 
ulation, thus the Joule heating to the accretion flow is 
neglected. Obviously the electrons will be heated to a 
higher temperature when this effect is considered. 
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